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Abstract 



The antideuteron and antihehum-3 production rates at high-energy heavy ion col- 

^ ' hsions are calculated in the framework of fusion mechanism. It is supposed, that 

0^ . p, n, d, ^He participating in the fusion are moving in the mean field of other fireball 

I constituents. It is demonstrated that at high energies, where many pious are present in 

^f) • the fireball, the number of produced d and ^He is determined by the balance between 

. created and desintegrated (mainly in collisions with pious) d and ^He. The explicit 

I formulae for coalescence parameters are presented and compared with the data. 



PACS numbers: 25.75.-q, 25.75.Dw, 12.38. Mh 



X : 1 Introduction 



In recent experiments on heavy ion collisions the production of light antinuclei — an- 
tideuterons and antihelium-3 — was measured [1],[2],[3],[4]. Once may expect, that this 
process proceeds at the intermediate stage of the evolution of the fireball, created in heavy 
ion collisions. Because of their small binding energy antinuclei are formed at the stage, 
when the hadrons are already formed, the density of hadronic matter is of order of normal 
nuclear density, and the particle collisions are still important. We will call this stage as the 
"dense gas" stage of fireball evolution. It is important to have a theoretical description of 
antinuclei production in heavy ions collision, because the dense gas stage is evidently after 
the stage, where quarks and gluons transform into hadrons, but on the other hand this stage 
precedes the last stage of the fireball evolution (sometimes called the thermal freeze-out), 
when hadronic spectra are formed. The theoretical calculations of antinuclei production and 
their comparison with experiment can shed some light on the dense gas stage of the fireball 
evolution and allows one to do one step back from the final stage, about which we have 
direct information from experiment. Particularely, as it will be shown, the volume of the 
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fireball and the liadron densities at the dense gas stage can be estimated from the data on 
antinuclei production. It is a common belief, that the antinuclei production proceeds as a 
fusion process: 

p + n ^ d (1) 

in case of antideuterons and 

p + p + n ^ ^TTe (2) 

in case of ^He. According to the dominant coalescence mechanism (1,2) it is convenient to 
characterize d and ^He production by the coalescence parameters 

B. = e/^(e/;^\\ p,=p,/2. (3) 



d^Pd \ d^Pp 

d^ f d'Np _ 



d^Pm V ^ d^pp 

where one assumes d^Nn/d^Pn — d^Np/d^pp. In what follows we consider only the central 
heavy ions collisions. 

The basic ideas of our approach are the following. We assume, that the coalescence 
mechanism (1,2) is the dominant one in the production of heavy nuclei. (For d production it 
will be shown by direct calculation, that the contribution of competing process p+n ^ d + n 
is small, see Appendix B.) The fusion reactions (1,2) cannot proceed if all particles are on 
mass shell. However, in the fireball at the dense gas of its evolution, p, n, d, ^He are not on 
mass shell, since they interact with surrounding matter. One may consider their movement 
as a propagation in the mean external complex field caused by the matter. The interaction 
with tliis field leads to appearance of the mass shifts and widths of all particles propagating in 
the medium (or width broadening for unstable ones) , analogous to refraction and attenuation 
indeces in the case of photon propagation. Another important ingredient of our approach 
is the balance of antinucleous production and desintegration rates. This balance is achieved 
because of large density of pions in the fireball and high rate of 71+ {d or ^He) collisions 
leading to antinucleous desintegration. The balance does not imply a statistical equilibrium, 
but rather a stationary process. The statistical or thermal equilibrium are not assumed in 
the calculation. 

In most of previous investigations the production of deuterons or light nuclei was con- 
sidered, but not the production of antinuclei. These processes have some common features, 
but also some differences. 

The production of light nuclei in heavy ion collisions nas been studied for many years (for 
an early review see [5]). The first calculation of the deuteron production rate was performed 
in [6] within simple model of optical potential in nuclei, resposiblc for the deuteron formation. 
In [7] it was proposed that the proton and neutron bind together if their relative momentum is 
less than some po, a phenomenological parameter to be determinad from experiments. Then 
B2 — ■ \ • ^Pq is proportional to the propability to find the proton and neutron inside 
the sphere of radius po in momentum space. A thermodynamical approach to the deuteron 
formation was applied in [8], where the deuterons were assumed to be in thermodynamical 
equlibrium with other nucleons in the fireball. 
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In [9] the sudden approximation of quantum mechanics was first apphed to the hght 
nuclei formation. It imphes short transition from the high-density stage \i) consisting of 
the protons and neutrons only, to the low-density stage |/), consisting of the deuterons as 
well. The amplitude of the deuteron formation is given by the overlap of the wave finctions 
Under simplifying assumptions that these states are free moving particles, uniformly 
distributed in a box of volume V and neglecting the deuteron size, one obtains the coalescense 
parameter: 

^ ^ ^ 3 (27r)3 



Ep 4 V 



where 3/4 is the ratio of the spin weights, the factor 2/Ep comes from invariant definition 
of the phase space element d^p/E. 

Most of investigations of the light cluster formation in heavy ion collisions represent 
extentions of the sudden approximation result for variously prepared initial state (sourse). 
Among them are different sourse geometries [10], [11], expanding systems [12], intranuclear 
cascade model [13], diagrammatic approach [14], relativistic quantum molecular dynamics 
model [15], Boltzmann-Uehling-Uhlenbeck transport model [16]. In most cases the results 
are based on computer simulations and analytical cxpcssions arc unavailable. 

A novel phenomenon is observed if one tries to go beyond the sudden aproximation. In 
[17] the deuteron formation was considered in the time-dependent perturbation theory. If 
the interaction, responsible for the deuteron formation, is switched on during some finite 
time T, the coalescense parameter gets multiplied on additional factor ~ (sr)^/^, where e is 
the deuteron binding energy. Exact value of r is not determined; as argued in [17], it should 
be of order of colliding nuclei radius, so (er)^/^ ~ 1. 

In all mentioned above approaches the interaction of the deuterons with the fireball 
enviroment was not accounted. As it will be shown below, this interaction, especially with 
pions, plays an essential role. As a consequence the coalescence parameter acquires the factor 
(s/r)^/^, where T is the deuteron in-medium width, which is calculated here. According to 
our calculation (e/r)^/^ _ ^/^q. 

Another transport model of the deuteron formation in nuclear collisions is available in 
literature [18, 19]. It was assumed, that the deuteron formation and desintegration goes 
mainly via triple nucleon reaction NNN ^ dN. At low (nonrelativistic) energies the re- 
action rates were calculated using impulse (Born) approximation with the help of Faddeev 
equation. Again, the thermodynamical equilibrium was assumed in the calculations. This 
approach also does not account the deuteron desintegration by the pions. 

It must be stressed that we apply our approach only to the process of antideuieion (or 
antihelium) formation in high-energy central heavy ion collisions. In this case, prior to 
antinuclei formation, the hadrons experience at least few collisions with each other and the 
concept of the hot dense gas seems reasonable to apply. Although the approach is valid for 
the deuteron (or helium) formation processes, the results should be used with care: a part of 
the outgoing deuterons may be just some fragments of primary colliding nuclei. Is is rather 
difficult to separate this part from the deuterons, synthesized from the hot hadronic gas, 
at least some experimental cuts should be applied. We will not consider this point here, 
restricting ourselves only by antinuclei processes. 

The material of the paper is presented as follows. In Section 2 the effective low energy 
Lagrangian for fusion processes (1,2) is constracted. The value of effective coupling constant 
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Figure 1: Elastic pn scattering through d formation 

is found by two ways: 1) by considering the elastic scattering amplitude p+n p+n (in case 
of d production); 2) by consideration of d, or "^He polarization operators, on the base of the 
field theory. With the help of effective Lagrangian the cross sections of fusion reactions are 
calculated. In Section 3 the transport equation for antinucleous production and propagation 
in the fireball is formulated. The formulae for the widths and mass shifts of particles, moving 
in the medium are presented and the formation rates of d and ^He are calculated in terms 
of effective p, n, d, ^He widths in the medium. In Section 4 the balance conditions for d and 
■^He production are formulated and the explicit formulae for coalescence parameters (3), (4) 
are presented. Section 5 is devoted to comparison with experimental data. The model of 
pion and nucleon densities at the dense gas stage of the fireball evolution is formulated. 
With the help of this model the values of the widths were calculated and it is demonstrated, 
that all assumptions, used in the calculation are fulfilled. The coalescence parameters are 
calculated numerically for experimental conditions and compared with the data. Section 6 
presents our conclusion. The details of the calculation are given in the Appendices. 



2 Effective low energy Lagrangians and fusion cross 
sections. 

At first let us consider the antideuteron production. Evidently, (7p+^_^j = ap-\-n-td- The low 
energy effective interaction Lagrangian has the form 

L = giJpipn^Pd , (5) 

where ipp, ipn and ipa are nonrelativistic p, n and d wave functions and spin effects are 
neglected. Let us temporarely suppose, that pnd interaction is point-like, i.e. neglect the 
nuclear force radius Tq. Following Landau [20], calculate the forward p + n — > p + n elastic 
scattering amplitude. In this approximation it corresponds to Fig.l diagram. According to 
standard rules the relativistic amplitude is given by: 

M = ^ fa 6 

s-ml Am(E + £) ^ ' 

where s = (pj, + Pri)^, m is nucleon mass, e — nip + m„ — = 2.225 MeV is the deuteron 
binding energy, E is the cm. energy, = rUp + m„ + E. 

On the other hand, the nonrelativistic s-wave amplitude / of elastic scattering according 
to Bethe-Peierls [21] theory is equal 

A;(ctg5o — k — ia^ 
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where a — y/me, k is the proton (or neutron) momentum in cm. frame. At the deuteron pole, 
s — > m^, k ia, the amphtudes (6) and (7) must be equal (keeping in mind the different 
normalizations M = STT^/sf). This requirement determines the value of the effective coupling 
constant [20]: 

/ = 128n m^/me (8) 

It must be stressed, that g'^ — > at £ — > 0, which means that the fusion rate goes to zero 
at £ ^ 0. Physically this is natural: the probability of production of large size antideuteron 
in heavy ion collision decreases with increasing of its size. It is an easy task to account for 
the finite value of nuclear forse radius tq, see [22] and Appendix A. As a result, the forward 
scattering amplitude /o at the pole k = ia gets an additional factor (1 — aro)~^ ~ 1.67. 
With account of this correction g'^ becomes equal: 

g'^ = 1287: m^/rne {1 - aro)~^ (9) 

The invariant cross section for the fusion process pn ^ d is found by common rules: 

E^f^ = r-r^=^=j\p.+Pn-P.), (10) 

where Ep,En and Ed are p, n and d total energies, 3/4 is the ratio of the spin weights. 
In fact, the cross-section (10) is valid only for low (nonrelativistic) cm energies. At high 
energies excited states are produced (for instance A^A resonance) which may decay into the 
antideuteron. Here wo do not consider these rather complicated (in medium) processes, 
assuming they are small enough. 

Let us now derive the expression for the effective coupling constant (in the limit Tq 0) 
by another method - by consideration of the deuteron Green function and polarization 
operator in field theory. This approach, of course, gives the same result as above in case 
of deuteron, but it will be useful in consideration of ^He production, where the scattering 
theory is rather complicated. 

For nonrenormalized deuteron Green function D[p'^) we have the Schwinger- Dyson equa- 
tion: 

[p' - ml, - U{p')]Dip') = 1, (11) 
where q is deuteron bare mass. After mass renormalization we get 

{P' - ml, - n(ml) - [n(/) - Tl{ml)]]D{p') = 1, 

<o + n(m2) = ml, (12) 

where is the physical deuteron mass. Then we perform the Green function renormaliza- 
tion -D(p^) = Z2Dj.en{p^) so that 

f d p^ — m^ 

The main assumption in this approach, used in derivation of the Bethe-Peierls equation (7), 
is: vertex corrections and corrections to nucleon propagators are neglected, d is assumed to 
be the bound state of p and n. The polarization operator is contributed only by the diagram 
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a) VJL-^ b) 

Figure 2: Deuteron and helium-3 self energies 

of Fig. 2a) and we can write U — QqII, where is the nonrenormahzed couphng constant, 
related to renormalized one g'^ as: 

g'o = z^V ■ (14) 

Representing U in (12) by dispersion relation we get 
The condition of (13) results in: 



Z2 + - ds- ^ = 1 16 



If deuteron is the bound state of proton and neutron, then the contribution of bare deuteron 
to Schwinger- Dyson equation vanishes and Z2 — 0. Putting Z2 = in (16) we get the 
equation for g^: 

Im n(s) can be easely calculated for Fig. 2a) diagram (nonrelativistic approximation can be 
used). As a result the same expression for g^ (8) as in the first method is obtained. 

Let us turn now to the study of ^He production. The effective nonrelativistic Lagrengian 
has the form 

L = Go^pi^p2^n^;^e (18) 

The derivation, performed above in the second approach, can be repeated and eq. (17) 
with substitution g"^ — > G^, rrid —>■ rriHe is found. The Imll{s) is described now by the 
diagram of Fig. 2b). There is, however, the essential difference in calculation of the Fig. 2b) 
diagram contribution to (17) in comparison with that of Fig. 2a): the integral in (17) is 
linearly ultraviolet divergent in nonrelativistic approximation. (This divergence still persists 
in relativistic calculations.) This circumstance physically corresponds to the well known fact, 
that in nuclear physics the 3-body problem with 5-function potential cannot be correctly 
formulated: it is necessary to have an additional information about the interaction at small 
distances. For our purposes it is sufficient to have a rather crude estimation of G^. So, we 
put ultraviolet cut-off in the integral in (17). The calculation gives: 

= 36^3 (47r)3^ (19) 

In numerical calculations A = 300 MeV will be taken. 
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3 Transport equation. Calculation of collisions inte- 
gral. 

As mentioned in the Introduction, we suppose, that antideuterons and antihehum-3 are 
formed at the dense gas stage of the fireball evolution, which follows after "chemical freeze- 
out" stage [23, 24] . Let us assume, that particle propagations at this stage can be described 
by means of classic transport (kinetic) equations. 

For definiteness, consider first the rf-production. We will use the notation qi{x,p), 
i = p,n,d,n, ... for the double densities in coordinate and momentum spaces and ni{x) = 
J qi{x,p) d^p for the coordinate densities. {qi{x,p) are Lorentz invariant.) Let us choose the 
cm. frame of colliding ions. The transport equation for qgix, p) has the form: 

md dqd{pa,x) ^ j _ dqg , ^ „ 

'e', dx, - -m 

= J d^Ppd^Pnqn{Pp)qn{Pn)uJpfi^d " Qd{Pd)Yl J d^Pi<iiiPi)^di^X (20) 

where — (1, vj)/y^l — v| is antideuteron 4- velocity, Up^^^d is the fusion reaction propa- 
bility, proportional to differential cross section 



- — e^e: d^ 

and similarly for the desintegration process di ^ X due to collisions of d with i-th constituent 
of the fireball {i = 'K,K,p.,n., etc.). The terms, where d appears in the momentum interval 
Pdi Pd + ^Pd ^ result of elastic collisions are neglected. The cross section ap^^d is defined 
by (10), the cross sections of d collisions with fireball constituents (cx^j, ...) can be found 
from the experimental data. Necessary applicabihty conditions of (20) is that particles 
wave lengths Aj = p^^ should be much less, than the mean distances d between fireball 
constituents, Aj ^ d. 

At the dense gas stage of the fireball evolution all particles inside the fireball should be 
considered as moving in the mean field of other fireball constituents. As a consequence, the 
masses are shifted in comparison with their vacuum values. Similarly, due to interaction 
with medium constituents, the widths Fj appear (or width broadening, if the particle has 
its proper widths). The mass shift Am{E) and width T{E) are expressed in terms of the 
forward scattering amplitude f{E) of the particle on i-th medium constituents (see[25],[26] 
and references therein) 

Am{E) = -27r^ Re fi{E) 

r{E) = 47r^/m/,(i?) = ^a,(^), (22) 
m m 

where E,p and m are particle energy, momentum and mass, Ui is the density of i-th con- 
stituent in medium. Eq's.(22) take place in the system, where the constituents are at rest. 
In the case of moving constituents the corresponding Lorentz boost must be performed. (By 
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definition Am and F are Lorentz invariant, for details see [27]). The relations (22) must be 
summed over all fireball constituents. 

The applicability conditions of eq. (22) are the following [25, 26]: 

i) . A < d 

ii) . 

iii) . The main part of the scattering proceeds at small angles, 6* -C 1. 

The conditions i) and iii) are well satisfied in the cases under consideration. The condition 
ii) is fulfilled not quite well and some corrections should be taken into account (see below). 

Therefore, p, n and d can be considered as Breit-Wigner resonances with varying masses 
distributed according to Breit-Wigner formula. In the process of the fireball expansion these 
Breit-Wigner resonances smoothly evolve to their stable counterparts. So, we substitute 
(10) in the first term in the right-hand side of eq. (20) and integrate over the masses m' of 
Breit-Wigner resonances. Using the notation I / Eg for this term we get: 

' - fdKdrr^drr^. ^^"^ ^^"^ ^^"^ 



Ed J P ^ d _ ^_)2 + r|/4 (mk - mnf + r|/4 [m';^ - mj)^ + f |/4 

""T^wJ^^ ?P-(PP-) S'ipp + p. - p,) S{E'^ + - E'^ (23) 

where E'^ — ^ m'^ -|- etc. Assume, that Fp n is much smaller than m, and the variation of 

(lp,h{p) on the interval F is also small enough, Tdq/dp ^ q. (In fact, Fp = F^j ^ 200 MeV, 
see below and Am ^ 30 MeV.) Then the distributions qp{pp) = qniPn) can be taken out 
from the integral sign at the values Pp = Pn = Pd/'^- With a good accuracy we can put 
Tp = F^ = F. But Fj in (23) generally is not equal to the antideuteron in-medium width 
Fj 2F. The reason is that pn system with d quantum numbers at high excitations does 
not necessarily evolve to d in the process of the fireball expansion, but can decay in other 
ways. One may expect Fj < Fj. We keep the ratio a = Fj/F^ as a free parameter in the 
calculation. As can be seen below, the results weakly depend in this ratio. The calculation 
of the integral / is straightforward and leads to 



37r^ , /F(l + a 



- 8 » V^'^'»' 

The contributions of direct processes p + p — > d + 7i~ , n + n —>■ d + 7i~^ and p + n ci + 7r° 
are small, they comprise not more than 20% alltogether and may be neglected within the 
accuracy of our calculation, see Appendix B. 

We restrict ourselves to consideration of low and intermediate transverse momentum 
Pp^ < 1 GeV. At higher p± many phenomena go into the play: radial and elliptic flow, steep 
decreasing of spectrum with increasing of p± etc. The consideration of these effects requires 
more refined treatment of the problem. 

Turn now to the calculation of the second term in the r.h.s. of the transport equation 
(20), corresponding to d desintegration rate. The probability uJdi^x is proportional to the 
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total cross section (T^j^_^x the same way as (21). So one may write the disintegration term 
in (20) as: 

J d^Piqi{Pi)i^di-*x = (25) 
where Fj has an obvious interpretation of the antideuteron width. It is equal to: 

rj = — V / ^ qM J {PdPiT - mlmf a^i^x (26) 
rud'^ J Ei V 

It is Lorentz invariant generalization of the width (22) and depends on the momentum of 
the antideuteron, moving in medium. The process of ellastic vrd scattering was neglected in 
transport equation. 

It must be mentioned, that eq. (25) corresponds to gaseous approximation, when the 
screening corrections are not accounted. In fact, the pion densities in the fireball at SPS 
or RHIC energies are such, that the account of screening corrections is necessary. Such 
calculation, based on Glauber theory, is presented in Appendix C. As a result, it was found 
that the Glauber correction at SPS or RHIC conditions reduces Tp approximately by 30%. 

4 The balance condition. The formula for coalescence 
parameter for d. 

Suppose, that the rate of antideuteron collisions with other constituents of the fireball re- 
sulting in antideuteron desintegration is much larger, than the rate of the fireball expansion. 
This happens in case of heavy nucleous collisions at high energies, when the fireball size at 
the dense gasc stage is large, because of large number of produced pions per miclcon. The 
large relative pionic density in the fireball leads to high desintegration rate of antidcutcrons 
in collisions with pions. (The collisions with nucleons are less essential in d desintegration.) 
In this case one may expect a balance: the antideuteron production rate (first term in the 
r.h.s. of (20)) is equal to its desintegration rate (second term in r.h.s. of (20)). The balance 
condition determines d-density: 



The momentum distribution d^N^/d^p^ entering in (3) is obtained from (25) by integration 
over the fireball volume 

d'NM ^ /^3^^_^^_^^ ^28) 



d^Pd 

Prom equations (3), (27), (28) and (9) we find the coalescence parameter 



Since the x-dependence of qp{pp, x) is not known, we replace (29) by: 



R^^ - 247r\ / s(l + a) 2 n^. 
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Figure 3: Fireball at the dense gas stage. The effective volume Veff is the half of the outer 
shell of thikness from which the antiprotons reach detector 



where V is the fireball volume, rip and n'^p are the mean and mean square p densitites in 
the fireball. (The coordinate dependence of vT is neglected). B^^ is Lorentz invariant, 
as it should be. The volume V in (30) can be understood as a mean value of the fireball 
volume at the stage, where on the one hand, the hadrons are already formed, i.e. the mean 
distances between them are larger than the confinement radius Rc ~ l/'m.p ~ (1/4) fm, but 
on the other hand, hadron interactions are still essential and the mean distances between 
the antinucleons are of order or larger than the deuteron size, so d could be formed in 
pn collisions. The antinucleon distributions np{r),nn(r) inside the fireball are essentially 
nonuniform: at the stage of the fireball evolution, preceding the dense gas stage, antinucleons 
strongly annihilated in the internal part of the fireball and in much less extent in its external 
layer of the thickness of order p(n) annihilation length For this reason ri^p/n^ may 

be essentially larger than 1. For the same reason the antinucleons and antideuterons from 
the backside of the fireball (with respect to the observer) are absorbed in the fireball and 
cannot reach detector, sec Fig. 3. Therefore, only one half of the fireball volume contributes 
to the number of registered p, n and d. The corresponding factor approximately equal to 2 
is accounted in (29). 

The width F can be calculated in one or another model of fireball evolution. However, 

since it enters in the expression for the coalescence parameter Bij^ (30) as a/T, it influences 
i?2^ not significantly. The same remark refers to not quite certain patametcr a. Therefore, 
the comparison with the data allows one to find the most essential parameter VifipY /ri^, 
which would make it possible to check various models of the fireball evolution. 



5 The model of dense gas stage of fireball evolution. 

We accept the following model for the dense gas stage of fireball evolution [26] . (A related 
model had been suggested long ago [28, 29]: it may be called Fermi-Pomeranhuk model). 
Neglect for a moment contributions of all particles except for nucleons and pions. Assume 
that at dense gas stage any participant - nucleon or pion occupies the volume vat or Wtt, 
respectively. Then 

- V - 1 + - V - 1 + ^^'^ 

where = l/fjv, (3 = v^^/vn, Qtt = A^^/A^at- For numerical estimations we take (5 = 
(rjr/rTv)^ = 0.55, where r^r = 0.66 fm and = 0.81 fm are pion and nucleon electric radii. 
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Experiment 


AT A /I /I fl 1 

iNA44 [ij 


o iAK \l\ 




\fsl A, GcV 


17.4 


130 


4.8 


Pp±, GeV 


0.55 


0.33 


0.17 


-E'pcm, GeV 


1.5 


1.05 


0.99 


# of "wounded" nucleons 


362 [31, 32] 


320 [33] 


350 [34] 




5.2 [31, 32] 


7 ± 1 [35] 


1.6 [34] 



Table 1: Experiments, where antinuclei were observed 



We take as a parameter varying in the interval from = 0.17fm ^ (normal nuclear 
density) to 0.30 fm"^ 

Check first the applicability of our approach. Choose = 0.24 fm^'^ and Qt^ = 5.2. 
(The latter value was found at NA44 experiment [30] at SPS.) We have un = 0.062 fm~^, 

= 0.325 fm^, n = + = 0.39 fm~^ and the mean distances between the fireball 
constituents is d — — 1.4 fm. Evidently, the conditions A^f = l/pjf <^ d and Rc <S 

d are well satisfied. Check now the balance conditions that the probability of deutcron 
desintegration exceeds the fireball expansion rate. The former is given by 2T{m/ Ep). In 
calculation of T we assume, that the hadronic spectra at the dense gas stage of fireball 
evolution differ not too much from the spectra at the final (freeze-out) stage. Since, as was 
explained above, rather crude estimation of F is sufficient for our purposes, we believe, that 
such approximations is satisfactory. Using (25) and the pion spectrum, presented in [32], F 
was found, F ^ 200 MeV (Glauber correction was accounted). The estimation for the fireball 
expansion is -u; ~ (l/5)fm"\ We have: 2F(m/^p) ^ 1.4 fm"^ > 0.2fm^\ {Ep ^ 1.5 GeV 
at NA44 experiment [1]). So, this condition is also fulfilled. Finally, check the condition ii) 
of Section 3, which is equivalent to the requirement Im f <^ d. Before taking into account 
the Glauber correction Im / 1 fm and this condition is not well satisfied. After accounting 
30% Glauber correction the situation improves, but not too much. For this reason the values 
of F, presented above, have a large (may be 50%) uncertainty. Since VT enters (30), the 
error in B^^ reduces to 25%. 

6 Comparison with data on d-production. 

The antideuteron production in heavy ion collisions was observed in NA44 experiment at 
CERN in Pb + Pb collisions at = 17 AA GcV [1], STAR experiment at RHIC in Au + Au 
collisions at ^/s = 130A GeV [2] and E864 experiment at AGS in Au + Pb collisions at 
^/s = 4.8/1 GeV [3]. The data of these experiments are presented in Table 1. 

Table 2 gives the values of coalescence parameters i?2^'', measured in these experiments 
and relevant centralities. It must be mentioned that Sg^^ rather strongly depends on the 
centrality. For example in NA44 experiment the results for 0-5% centrality are about 1.5 
times lower. Also, B2^^^ has significant Ppt dependence: such dependence was observed at 
PHENIX experiment at RHIC (we do not analyse PHENIX data here, since only preliminary 
results were pubhshed up to date [4] ) . 

Table 2 presents also the results of our calculation for the values of n%, introduced in 
(31): n% — 0.17, 0.24, 030 fm~^. In the calculation of mean volume V at the dense gas 
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Experiment 


IN A44 


o lAK 


TT^/^O A 

hvo84 


Centrality, % 


n in 


n 18 
io 


n in 




Lie V 


A A _L 1 O 

4.4 ± V.o 


4.5 ± 0.3 ± 1.0 


41 ± 29 ± 23 




r, MeV 


148 


159 


109 


n% = 0.17 fm~^ 


V, lO^fm^ 


8.51 


10.95 


4.26 
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3.17 
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11.2 




r, MeV 
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206 


142 


n% = 0.24 fm"^ 


V, W' fm^ 


O.03 


7.7b 


O AO 

3.02 




Bf, lO-^GeV^ 


4.00 


4.22 


13.9 




r, MeV 


214 


232 


187 


n% = 0.30 fm"^ 


V, 10=^ fm-'^ 


4.26 


6.21 


2.42 




Bf, lO-^GeV^ 


5.28 


4.96 


15.0 



Table 2: Antinucleon width, fireball volume and antideuteron coalescense parameters 

stage the number of "wounded" nucleons and the number of pions were taken from Table 1, 
the corrections 15%, 20% and 10% were accounted for other particles, except for nucleons 
and pions, in cases of NA44, STAR and E684 experiments respectively. The antiproton 
annihilation length was estimated as Ip^ann ~ 3fm in case of NA44 and STAR experiments 
and Ip^ann ~ 1-5 fm in case of E684. In all cases we can put with sufficient accuracy n^/n^ ^ 2 
and a — 1/2. 

The calculated values of the coalescence parameters Bl'^ are presented in Table 2. It must 
be mentioned, that in case of E684 experiment the validity conditions of our approach are on 
the edge of their applicability. So, the theoretical expectations for B2 in this case are valid 
only by the order of magnitude. In two other cases the agreement of the theory with the 
experiment is quite satisfactory. Let us remark, that the main uncertainty arises from the 
effective width F which, probably, is known with an accuracy of order 50%. Taking in mind 
all other possible uncertainties, we believe that the accuracy of theoretical predictions for 
B2 is also about 50%. The inspection of the Table 2 shows, that the total hadronic density 
(pions+nucleons) at the dense gas stage of fireball expansion, where the antideuterons are 
formed, is of order un + JT-tt ~ 0.4 fm~^. Much lower or much higher densities would lead to 
contradiction with experiment. 

Using the pion and nuclcon densities and spectra chosen above, we can estimate the 
typical energy density at the dense gas stage: e ~ 0.3-^0.5 GeV/fm^. As the colliding ions are 
initially strongly Lorentz-compressed and the produced fireball expands in the longitudinal 
direction with an almost velocity of light, it takes approximately i? lOfm to reach the 
proposed volume. It is instructive to put these values on the Fig. 27 from [36], where the 
density of energy/time evolution of the fireball is represented. We see that the dense gas 
stage begins just after the moment when the total hadronization happened and this is a 
good check for self-consistency of our assumptions. 
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7 Production of ^He. 



The calculation of ^He production proceeds along the same lines as d. The only difference 
is that now the collision integral / - the first term in the r.h.s. of eq. (20) corresponds 
to the formation of ^He in collision of three antinucleons: pi,p2 and n and the effective 
coupling constant is given by (19). Instead of Breit-Wigner off-shell mass distributions for 
antiparticles pi,p2,n and ^He we now take the Gaussian ones: 

Mm',) = 2 ^32) 

The reason is that the Breit-Wigner distributions do not provide the necessary convergence 
of collision integral. The collision integral takes the form: 

/ TtG^ f f TJ j / £ r /\\ f ^^P^ ^^P^ ^^P 



dm'Ji{m'i)j J ^ ^ Qpipi) qpip2) QniPn) 



Ehp G4:Efjp I < -^-^ J I h/^ 

ne tie j j=i 2,3 Pi P2 Pn 

X 5\p, +P2+Pn- PHe) S{E'^, + E'^^ + E'^^ - E'^J (33) 

(The distribution of ^He was taken as a 5-function, what is equivalent to a = I" He/ (3r) <^ 1). 
The integral is equal to: 

G^r^'ip) (34) 



96^6 

Applying the balance condition and (19) we find the coalescence parameter: 

1 r 1 

Bs = 967r^ Ae-t^,^ (35) 

The accuracy of calculation of is lower, that in case of B2 since additional uncertainties 
appear: strong dependence on ultraviolet cut-off A etc. So we may pretend only on the 
estimation of B3, correct up to order of magnitude (in the best case up to factor of 2). The 
coalescence parameter for '''He production in Au + Au collisions at ^/s = 130A GeV was 
measured by STAR Collaboration [2]. In this experiment the centrality was up to 18%, ^He 
transverse momentum 1.0 < p^^^ < 5.0 GeV and rapidity \y\ < 0.8. These limititations 
corresponds to the average antiproton energy Ep f« 1.4 GeV. STAR found: 

Bl'^lj- = (2.1 ± 0.6 ± 0.6) X 10"^ GeV^ (36) 

The theoretical value according to (35) at n% = 0.24 fm"'^, rfi/v? = 3, and A = 300 MeV is: 

B^ jj^ = 3.3 X 10-^GeV^ (37) 

The agreement with experiment is good, despite of many theoretical uncertainties. It demon- 
strates the vahdity of basic ideas of theoretical approach. 
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8 Summary and conclusion. 



The coalescence parameters for production of antideutrons and antihelium-3 in heavy ion 
coUisions were calculated. The obtained results are based on three assumptions: i) the main 
mechanism of light antinucleous production is coalescence (fusion) mechanism — eq. (1,2); 
ii) all particles, participating in fusion process are moving in the mean field of other fire- 
ball constituents; iii) the number of produced antinucleous is determined by the balance 
conditions: the equality of produced and desintegrated — mainly by pions — antinuclei. 
The production of antinucleous proceeds at the dense gas stage of the fireball evolution, 
when the hadrons are already formed, but their interaction is still important. Statistical 
or thermal equilibrium are not supposed at the dense gas stage. In fact, the final results 
depend on one parameter — the volume of the fireball at this stage (or, equivalently, on 
the hadron densities.) Good agreement with experimental data for coalescence parameters 
was obtained for experiments at CERN, RHIC and AGS for the values of n%, defined by 
eq. (31), n% ~ 0.17 -r 0.30 fm~^, close to normal nucleous density. Much lower, or much 
higher values of n% lead the to values of coalescence parameters, incompatibe with the data. 
It must be stressed, that the same values of n% well describe the coalescence parameter 
for ^He production, which demonstrates the effectiveness of the method. The values of the 
fireball volume V are about 2 times larger, than those found in [23] at the so-called chemical 
freeze-out stage, about 2 times smaller than at thermal freeze-out [35, 37] and in agreement 
with limitations on the volumes found in [24]. The same method can be applied to the 
production of deutcrons and ^He nucleous. But in these cases, in order to avoid the back- 
ground, consisting of d or ^He falling apart from colliding nuclei, it is necessary to select the 
events with rapidity close to zero. One may expect, that for such events the expressions for 
coalescence parameters, found above are valid. The further investigations of this problem — 
both theoretical and experimental are very desirable, since they can shed light also on the 
most interesting stage — the stage of hadron formation. 
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Appendix A: Account of the finite nuclear force radius 

According to Bethe the general form of pn scattering phase 5 at low energies is 



Here k is the proton (or neutron) momentum in cm. frame, tq is defined as a radius of 
nuclear forces: at r > rg the pn potential Vpn{r) may be neglected. (A.l) is an expansion in 
terms of krQ, higher order terms, ~ k^r^ are omitted. For the scattering amplitude (7) we 




(A.l) 
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get from (A.l) 

J. 1 i 1 (A 2) 

A;(ctg5o — k — ia 1 + i{k + ia)ro/2 

The value of effective coupling constant is determined by the residue of the amplitude / at 
the pole k — > ia. It is easy to see, that in comparison with the case ro = 0, the value of the 
residue changes by the factor (1 — aro)~^. Higher order corrections ~ k'^rg are negligible. 

There arc several different definitions of the nuclear force radius Tq. On the other hand, 
all rp-corrections to the pnrf-coupling can be expressed in terms of the normalization constant 
As of the deuteron radial wave function R = Ase~^^'^ at large r. Indeed, it determines the 
residue of the amphtude (7) at the pole E — —e: 

In case of zero radius of nuclear force As = (4me)^/^ = 0.68 fm^^^. Analisys of experimental 
data on elastic pd scattering [38] as well as various potential models give the value around 
= 0.88 fm^^^. Comparison with the amplitude (7) gives the following coupling value: 

= 647rmA| (A.4) 

It exceeds the zero radius result (9) by 1.67 times. It corresponds to Tq = 1.7 fm. 



Appendix B: Estimation of direct process AT AT diz 

The main processes of the deutron formation in vacuum is NN — > dir. In case of exact 
isospin symmetry the total cross sections of each chanell are related as 

similarly for antinucleons. The process pp dn^, as well as inverse one, has been accurately 
measured in many experiments. The cross section is shown in Fig 4. It has a resonanse peak 
near ^/s = 2.17 GeV, the mass of the NA system and decreaces steeply above s ~ 2.5 GeV. 

In this appendix we estimate the contribution of these processes to the deutron formation 
rate. It is convenient to use invariant definition of the cross-section in relativistic calculations. 
We assume qp{p) = (Inip)- The processes with identical particles in initial state pp —>■ dir'^, 
nn d7r~ should be taken with weight 1/2 to avoid the double counting from the same 
regions of the phase space. So the total deuteron production rate can be written as: 



"^d u9qd{Pd) 3 /■ 3 3 3 / \ / \ hirect 

~ir^d g^f, ^ 2 / ^ P"^^ P''Wv%\P2)^pp^d^+ = ' ^^-^^ 

where uj is the reaction probability: 

_ V ^^pp^dn+ (-r, Q\ 

E1E2 d^pdd^p^ 



Eh dx^" 2 
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2 2.1 2.2 2.3 2.4 

Figure 4: Total cross section app^d-K+ in mb versus ^/s in GeV. The experimental data are 
taken from [39]. The line represents the fit of [40]. 

The value Idirect should be compared with the collision integral / of the fusion process, 
obtained in (23,24). Integrating by the momenta p^r) one obtains: 

hirect = \ j (ip(P^) (ip(P^) \l (PiP^y -^^e / ^'^^^'^^^ (B.4) 

We compute this integral in the center of mass frame. Let (£, k) and (f^^, k^) be the cm. en- 
ergies and momenta of the first proton and the deuteron respectively. Then the lab momenta 
are given by the Lorentz boost, v is the center of mass velocity: 



vfvk) 

k + (7 - 1)^^ + 7v^ , = 7(£ + vk) 



Pi " ■ V/ 2 

V 

P2 = -k-(7-l)^^^ + 7v£, £;2 = 7(^-vk) (B.5) 

where 7 = (1 — v^)'^!"^ . Now in the integral (B.4) we come from the variables pi, p2 to v, k 
according to: 



El E2 

Notice, that Ed^cr/d^p is invariant cross-section, so E^d^a/d^pd — Sdd^a/d^kd ■ Then the 
integral (B.4) can be written as: 

h^rect = \ -^^j d^vd^kq^ q,{p,) ^^£^k£d^^^^ (B.6) 

And finally, it could be convenient instead of the velocity v to substitute the deuteron cm 
momentum k^ according to: 

2w pd - kd 
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The velocity integration measure becomes: 

4,3 8 ^3 f Ed + Sd Yd^kd 

Then the integral (B.6) can be written as follows: 



' direct 



2 J d^kd {mj + EdSd + Pd^df 



To evaluate it further, one needs to know the antiproton distribution qp{p) in the fireball. 
For the first approximation, we could take qp{p) ~ QpiPd/'^), aw we did in (24). This however 
would be an overestimation of the production rate, since it gets the main contribution from 
the resonace area (see Fig. 4), but not from the threshold. More carefuU estimation can be 
obtained by taking the Boltzman distribution for the antiprotons: 

qp(p) ~ e-^/^ (B.IO) 

Then one can perform the integration over the angles and kd- The result is: 

hired = qliPd/^)—^ / ds ■ s{s - Ami) app^d7T+ 

( Edjs -ml- ml) \ sh{^pdkd/{Tml)) 
" I ^d ) V-sPdkd/iTml) (^-^^^ 

where app^dTr+ is the total cross section, k^—[s — {md — m.j^)'^][s — {md + rriT^Y]/ (4s). 

The integral (B.ll) reaches maximum for the dcuteron at rest in the frame, where the 
system is described by Boltzman distribution (B.IO). So for numerical estimation we take 
Pd — ^ and the temperature T of order of the antiproton inverse slope parameter T = 
300 MeV. We intergate (B.ll) up to ^/s — 2.5 GeV. Above this value other channels appear, 
where experimental data are rare. However the total cross section pp — > dX at high ^/s is 
neghgible, so this will not be an underestimation. The result is 

Idirect = 7.5 GeV^ X g|(pj/2) (B.12) 

It should be compared with the fusion rate integral (24). According to our estimations 
I = 40 GeV^ X g|(pj/2). It confirms our assumption that the contribution of direct process 
of antideuteron formation NN — > dir is small in the heavy ion collisions. 



Appendix C: Particle width in dense medium 

The relation (22) for Tp corresponds to the case of nucleon moving through the medium 
of pions at the dense gas stage of fireball evolution. It turns out that Tp weakly depends 
on the pion spectrum details, so with sufficient accuracy one can use the observed pion 
spectrum which is formed on the freeze-out stage. To obtain numerical results we exploit 
eq. (26) for pion-antiproton collisions, where the pion distribution is extracted from [32]. 
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It is usually parametrized in terms of transverse mass m^J: = •y'^(pj)^~+mf and rapidity 
y^\H{E + p\\)/{E-p\\)]: 

27rE^q^{p) = ^ \ = Cexp{-{m^ - m^)/T^ - y^jh^}, (C.l) 

where Tj_ = 110 MeV, h = 1.9 and C ~ 38.6 GeV^^ is the normalization constant, determined 
from the requirement J qT^{p)d^p — Ut^. Performing the calculation we find that the numerical 
value of Tp is about 290 MeV at the pion density = 0.32 fm~^ and antiproton energy 
Bp = 1.5 GeV. The effect of averaging over pion momentum is practically negligible as the 
obtained value of Tp differs by less than 5% from its value calculated within the assumption 
that all pions are at rest. 

As it was pointed in the text, Eq. (22) imphes that the scattering amplitude of the 
incident particles on the medium constituents is much less then the mean distance between 
them and the scattering can be treated as the sum of independent processes. However, this 
condition is not fulfilled quite well in the dense matter. In this case na has the meaning 
of the total cross section on all the constituents of fireball (divided by its volume) which is 
not simply reduced to the sum of individual cross sections. This lack of additivity is well 
known from the nucleon-deuteron scattering at high energies [41]. Within the accuracy of 
our approach it is sufficient to take only the first term in density expansion. Consider the 
scattering amplitude at small angles when the incident nucleon interacts with the target 
consisting from two pions confined in a cube with side a. From the calculation of Tp it 
follows that the main contribution to the integral for Fp arises from the pions with small 
momenta, so one can assume that the antiproton momentum is much larger than the pion 
momentum. Then the antiproton scattering amplitude can be written in the framework of 
method developed by Glauber [41]: 

^^'^^ ^^ij d'pe-^''^^^^[eM'^tS,ip - fi) + 2t6,ip - f^)} - 1], (C.2) 

where 5i is the phase shift induced by the scattering of antiproton on i-th pion and is the 
transverse coordinate of i-th pion. 
Using the relation 

2; 



exp{2z5(r) - 1} = ^ / ^^^^"^"^''(1^ ^^'^^ 



and identity e^^"*"^^ — 1 = (e^ — 1) + (e^ — 1) + (e^ — l)(e^ — 1), we get the following equation 
for the forward scattering amplitude: 

/(O) = /i(0) +/2(0) + ^Jim2i0). (C.4) 

To obtain quantitative result we assume that /i = /2 and Re fi = (actually. Re f / Im f ~ 
20%) and note that n — 2/a^ which gives the first correction to na: 



n^(7 — > n7r(7 



4a2 V 2 



(C.5) 
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The factor in square brackets represents the screening effect and also appears in Fp. At the 
pion density n^r = 0.32 fm~^ and a^p ^ 30 fm^ the screening effect reduces the value of Tp 
to 190 MeV. (It is worth to note that the screening correction (C.5) to is larger than the 
same one to Tp as cr^j > cr^p which is a direct indication that Fj < 2Fp.) Taking different 
values of the pion density one can find the other values of Fp presented in Table 2. 
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